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Within a two-band tight binding model, we investigate the dynamics of electrons with Markofnan 



(N 



dephasing under the influence of static electric fields. With the help of both numerical and analytic 
calculations we find that the dephasing ultimately takes electrons which are initially located in one 
miniband to equal population of the two minibands, instead of undergoing persistent Rabi flop, as 
they do in the absence of scattering. Miniband localization is wholly destroyed by the intervention of 
dephasing. We also obtain the effective decay time for the approach to equal band populations under 
*^ conditions of small interband communication and in the long-time limit, through a perturbative 

r/2 ' calculation. The decay rate shows characteristic sharp peaks at values of the parameters which give 

D Zener resonances. 
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I. INTRODUCTION 



Recently there has been intensive investigation 1 7 of the influence of external electric fields on semiconductor 
superlattices. The availability of high-quality semiconductor superlattices, the possible practical application of new 
phenomena, and the surge of new interesting theoretical predictions have-attracted both theoretical and experimental 
attention to this active field. An important example is Bloch oscillationa (BO), one of_the earliest predictions of the 
band theory of solids. Its observation in standard crystalline solids has been prevented^ by excessive scattering under 
all practical conditions. But in semiconductor superlattices BOs have not only been observed experimentallyo but have 
been further exploited to make an extremely fast emitter of electromagnetic radiation — the "Bloch oscillator"!! Many 
other new phenomena have been predicted theore. ;ically and/or observed experimentally, including negative differential 
. conductivity™!, absolute negative conduct anceE , inverse Bloch— oscillation (strong THz=photocurrent resonance at 
' the Bloch, frequency)! 7 ], dynamic localizationErEJ, band collapsetJiilJ, band suppressionEZl, fractional Wannier-Stark 
0*N Laddersta, and multi-photon absorptiontH 

From a theoretical point of view the simplest models involve only a single miniband. But these miss all interesting 
interband phenomena, which are captured most simply in a two band tight binding model first introducedcS by 
G Fukuyama, Bari and Fogedby. Thisjppdel contains the essential physics of interband transitions, and can be compared, 
in practice with realistic situationsEj~c3 where only a single pair of bands is important. Rotvig, Jauho and SmitbO'CJ 
have studied coherent transport of one-dimensional semiconductor superlattices within this model. They found that 
coherent oscillation between the minibands can occur at special values of the applied field, where there ar£_.avoided 
crossings of the two interpenetrating Wannier-Stark Ladders (WSL) arising essentially from different bandso. These 
are the so-called Zener resonances. For electric field values between Zener resonances there are stable plateaus, where 
the relative population of the two bands is only weakly dependent on time. In each plateau there are small amplitude 
Bloch oscillations. The number of these, N, corresponds to the Nth Zener resonanceEa (labelled with decreasing field, 



so that the Nth resonance arises from aligning adjacent WSL levels in quantum wells N lattice spacings apart). 

In this paper, using the same model, we study the effects of scattering from lattice imperfections on the dynamics of 
electrons, using a stochastic Liouville equation for the density matrix, with Markoffian dephasing. By means of both 
analytic and numerical calculations we analyse the effects of scattering and miniband structure parameters_OB-,the 
time evolution of p_, the difference in electron population of the two minibands. In contrast to previous workc-iEj, in 
which the single quasimomentum component p~{k, t) was studied, we focus on the trace quantity p~(t) = p~(k, t). 
There are two reasons for this choice. The first is that the temporal evolution of the sum over all quasimomenta is of 
direct physical interest, describing the dynamics of the full set of electrons. The individual p_(k,t) depend not only 
on time but on quasimomentum k, so we should integrate p-(k,t) over the whole Brillouin zone. The other is that, 



with scattering, k is no longer a good quantum number — even in a suitable transverse gauge (see Ref. 22 and Eqs. 
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(17) - (19) below), where k is a good quantum number in the perfect lattice. We find that with no dissipation p-(t) 
periodically oscillates between 1 and — p m ("Rabi flop" between the initial condition of all electrons in one band and 
a fixed maximum exchange of a fraction of that population with the other band: < p m < 1) at avoided crossings, 
whereas when dephasing is introduced, the electrons initially assumed to be located in one miniband will ultimately 
be distributed equally between the two bands (band localization will be completely destroyed). 

In section II, we present the model and derive the equation of motion in k space. Section III provides numerical 
analysis of the effects of the single relaxation rate and of the miniband parameters on the Zener resonances. Section 
IV contains the perturbative calculation of the density matrix and its limiting long-time behavior. We make some 
concluding remarks in section V. 

II. MODEL 

We consider the standard tight-binding model of a two-band system in a static electric field E. The model 
Hamiltoniarc3 can be written as 

H = } (A a + nw B )a]ja„ + (A b + rvj) B )b\j3 n 

-(W a /4)(ai +1 a n + h.c.) + (W b /4)(bl +1 b n + h.c.) 

+eER(aib n + bla n ) . (1) 

Here the subscripts label the lattice sites and the lower and upper minibands are designated by symbols a and 6, 
respectively. We have introduced the notation lob = eEd for the Bloch frequency, which will appear often below. The 
first two terms describe the site energies of the Wannier states in the presence of the electric field, and W a ,b are the 
widths of the isolated (E — 0) minibands induced by nearest neighbor hopping: e° ,b (fc) = A a ,b =p (W a .b/2) cos(fcd), 
where d is the lattice constant. The last term is the on-site electric dipolc coupling between minibands; eR is the 
corresponding dipole moment. This Hamiltonian does neglect Coulomb interactions and electric.jdkiLole elements 
between Wannier states on different sites, but it contains the essential physics for the problcmE90c3n3. Note that 
the hopping parameters W a ^b are written here with opposite signs, so that with both parameters positive the band 
structure at E = is of the standard nearly free electron character, with direct band gaps at the zone boundary. But 
the calculation to foliw is valid for arbitrary signs of the parameters. 

It is easily shownEd'Kj that the exact spectrum of H is two interpenetrating Wannier-Stark Ladders. But what do 
the corresponding states represent in terms of the occupation of the original bands as a function of time, and what is 
the influence of scattering? For vanishing dipole matrix element between bands, R = 0, there is no interband mixing. 
Each of the two bands gives rise to a single WSL. Clearly, when the electric field amplitude is such that the ladders 
become degenerate, even small values of R lead to strong interband mixing. The crossing of the ladders is "avoided" 
by any finite R, and the behavior at those avoided crossings (the Zener resonances) is of particular importance and 
interest. 

We start by defining the density matrix in the representation of the two bands, 

P(t) = E P^ntM, (2) 
ijmn 

where i,j = 1 or 2 are band indices: £iq (Cm) an d Cm (Cm) designate aj^ (a m ) and 6 m (b m ), respectively. Within a 
constant relaxation rate approximationc£l, the density matrix p(t) satisfies the following stochastic Liouville equation 
(SLE) (we set % = 1 throughout this paper), 

i^ = [H,p(t)]-iTp(t). (3) 
Here Tp describes the relaxation of the off-diagonal elements of p through dephasing: 

Vp=Y d {l-8 ij 8 mn )(a ij i^i). (4) 

ijmn 

The utility of this simplest form of the SLE has been discussed by Kenkre and collaborators (see Ref . |2^ and references 
therein). The parameters Qfy measure the loss of phase coherence between sites, or the scattering lifetime of band states 
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labeled by quasimomentum. The relaxation to zero of the lattice site off-diagonal elementii-pf p is then effectively 
an infinite temperature approximation. The finite temperature corrections are essentialESEH to understanding the 
current response to an electric field, a response which vanishes in the absence of a difference in thermal equilibrium 
population between states of different energy, but those corrections should not play a major role in the phenomena 
studied here. We note that in this limit there is no need to distinguishEZl between the relaxation rates of energy and 
momentum. 

Since we are interested in the dynamics of occupation of various band states, it is convenient to work in a wave 
vector basis, by Fourier transforming the density matrix. In general, since p mn is not translationally invariant (a 
function only of m — n), we have a full set p l £ = ^2 mn p mn (0 ex P[ — ikm + iqn] of Fourier components. But we will 

be interested in the wave vector diagonal band occupation numbers p]^ k (t) = p tJ (k 1 t). These evolve according to the 
corresponding Fourier transform of the SLE: 



i^p 11 (k, t) = tu B ^-p 11 (k, t) - eER[p 12 (k, t) ~ p 21 (k, t)] 



[ 2lT dk' 

-ia nP ll (k,t) + ia 11 — p n (k',t), 
Jo 27T 



(5) 



i—p 22 (k,t) = iujb q^P 22 (k,t) + eER[p 12 (k,t)- p 21 (k, t)} 

f 2v dk' 

-ia 22 p 22 (k, t) + ia 22 / -^—p 22 (k\t) , 
Jo 27T 



(6) 



ig- t P 12 (k,t) = iu B —p 12 (k,t) - eER[p n (k,t) - p 22 (k,t)] 
+ (A - W cos k) p 12 (k, t) - ia 12 p 12 [k, t) , 



(7) 



i-p 21 (k,t) = iuj B —p 21 (k,t) + eER[p n (k, t) - p 22 (k, t)} 



-[A-W cos kj p 21 (k, t) - ia 2 ip 21 (k, t) , 

where we have adopted the notation A = A a — A& and W = (W a + Wb)/2. By introducing 

p+(k,t) = p n (k,t) + p 22 (k,t) , 

p^(k,t)=p 11 (k,t)-p 22 (k,t) , 

p+-(k,t) = p 12 (k,t) + p 21 (k,t) , 

p- + (k,t) = ilp 21 (k,t)-p 12 (k,t)} , 



(8) 



(9) 
(10) 

(11) 
(12) 



and, for simplicity, taking an = a 22 = Oi\, ot\ 2 = a 2 \ = a 2 to reduce the number of parameters in the theory, we 
obtain 



d . . d . . 

— p+(k, t) - uj B —p + (k, t) = -ai 



p+(k,t) - 



2 " dk[ 
2vr 



p+(k',t) 



(13) 



—p-(k,t) - wb— p-(k,t) = -2eERp-+(k,t) -ai 



p-(k,t)~ I ^- p-(k',t) 

Z7T 



(14) 



d d ( \ 

— p+_(M) - vb-qtP+- (M) = - W cos k J p- + (k,t) - a 2 p+-{k,t) , 
J^p_+(M) - u B ^rP-+(k ) t) = -(A -Wcosk)p + _(k,t) + 2eERp_(k,t) - a 2 p-+(k,t) 



(15) 
(16) 
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The equation for p + (k,t) describes particle conservation. It is decoupled from the others, and we will ignore it 

in the following discussion. The equations for p-(k,d\, p^ (k,t), and p \-(k,t) can be reduced to the following 

ordinary differential equations in an accelerated basisEa, k(t) = k — Lost or, equivalently, in the transverse or vector 
gauge discussed above and in Ref. ^4 

^_X{k,t) = -2eER Z{k,t) - a x [X{k,t) - X Q (t)\, (17) 
—Y (k, t) = [A - W cos(fc - u B t)]Z(k, t) - a 2 Y(k, t), (18) 

—Z(k,t) = -[A - W cos(k — cjBt)]Y(k, t) + 2eER X(k,t) - a 2 Z(k,t) . (19) 

HereX(fc,i) = p_(k-u B t,t), Y(k,t) = p + -(k-co B t,t), Z(k,t) = p- + (k-co B t,t), mid X {t) = X(k, t)dk/2ir = 
P—(t), . In the absence of any scheme for solving these coupled equations analytically, we turn to numerical solutions. 

III. ZENER RESONANCE WITH DISSIPATION: NUMERICAL RESULTS 

In this section we focus on the influence of dissipation on p_(t), the difference in population between the two 
minibands. For further simplicity the relaxation rates a% and a 2 have been set equal, a% = a 2 = a, in the following. 

We emphasize again that due to scattering, the time dependences of X(k,t),Y(k,t) and Z{k,t) as described in 
Eqs. (17) - (19) depend on the sum of these. quantities over all values of the wave vector, whereas in the absence 
of scattering, each value of k is independentO'LJ. As we have discussed abc#e, we are particularly interested in the 
dynamics of the difference in electron population between the two minibandsEi as. described by p_ (t) . In the perfect 
lattice some one-component diagonal terms X(k,t) undergo Rabi oscillationE3~t3 — i.e., X(k,t) oscillates between 
+1 and —1 with a period approximately determined by the value given in Ref. 23, if we choose the electric field at a 
value corresponding to one of the avoided crossings of the two interpenetrating Wannier-Stark Ladders. But further 
analysis shows that for some other wave vectors X(k,t) will not undergo complete Rabi oscillations (will not reach 
the full extreme values of ±1), even if we choose the electric field to be at the avoided crossings. This is depicted 
in the lower panels of both Figs. 1 and 2, which differ only by the choice of system parameters, listed in the figure 
captions. The upper panel in each corresponds to k = and the lower panel to k = n. It can be clearly seen that 
X{k = 7r, t) falls short of the minimum value of -1 in both cases. But, for the parameters of Fig. 2, the minimum is 
more closely approached than it is in Fig. 1, although in Fig. 2 interband communication has been decreased. This 
can be understood by making the coordinate transformation t = t' + k/ujB in Eqs. (17) - (19) ( with no relaxation: 
a = 0). Then all quantities depend only on the single parameter t'\ X(k,t) = X(t'),Y(k,t) — Y(t'), Z(k,t) = Z(t'). 
The initial conditions are set at a time t and therefore are still k dependent in terms of t'. But when E is very large, 
we have t w t', so that X(k,t), Y(k,t), and Z(k,t) will be independent of k. Then all X(k,t) will approximately 
undergo full Rabi oscillations. 

In general this is the dominant reason that p~(t), which represents a sum over all wave vectors, oscillates between 
certain negative values (but never reaches —1 ) and +1, no matter how we choose the electric field. We demonstrate 
this in Fig. 3a, where we use the same parameters as those of Fig. 1. In this figure, Zener tunneling and the more 
rapid BO can both still be clearly identified, but some oscillations in the stable plateaus have been smeared out. 
In particular, when p-{t) reaches +1, only a single sharp cusp has survived. The smearing is due to interference 
between the BO's for different values of k, which appear in slightly different positions, as can be seen by comparing 
the upper and lower panels of Figs. 1 and 2. This interference is also a factor (but not the main one) in suppressing 
the minimum of p—(t) from —1 at Zener resonances. 

With the introduction of dissipation the dynamics change dramatically (even in the regime where the scattering 
rate iSpless than the Bloch frequency, cx/lob < 1, so that some BO's remain well defined; we will treat only this 
regimeEj). We demonstrate this in Fig. 3b, where we use the same parameters as in Fig. 3a except that a/ajs = 0.02 
instead of zero. In this figure, we see that p_ (t) still oscillates between negative and positive values, but the oscillation 
amplitude decreases strikingly with increasing time. Plateaus can still be identified in several earlier peaks of the 
envelopes, but the BO's are increasingly suppressed in successive plateaus. If we choose the ratio a/ws larger ( e.g., 
a/u>B = 0.3 in Fig. 3c, where we use the same other parameters as in Fig. 3a), a few BO's can be seen but p-(t) 
never becomes negative. This demonstrates that observation of Rabi oscillations is in general more difficult than that 
of BO. 
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In the absence of scattering we can localize the electrons in one of the minibands simply by choosing the electric 
field amplitude to be well away from the avoided crossings. This can be seen through numerical calculation of P-(i). 
However, the introduction of dephasing destroys this picture completely. For any value of the electric field the electron 
density will be distributed equally between the two minibands at long times. This can be seen clearly by comparing 
Figs. 4a and 4b, where we use the same superlattice parameters as in Fig. 2, and the BO frequency u>b has been 
changed from the value at an avoided crossing of 10.2 meV to 9 meV. In Fig. 4a there is no dissipation, while in Fig. 
4b, a/ujB = 0.35. 

IV. PERTURBATION THEORY AND LONG TIME TRANSITIONS: EXACT RESULTS 

Under certain conditions pcrturbative calculations can give some useful results. Setting A(k,t) = Y{k,t) + 

iZ(k,t), B(k,t) = Y(k,t) — iZ(k,t),/j, = 2eER, and recalling a\ = a 2 = a, we can rewrite the Eqs. (18) and 
(19) as 

-^A(M) = -»[A - Wcos(k - Lj B t)]A(k,t) - aA(k,t) + ifj,X(k,i), (20) 

j t B{k, t)=i[A-W cos(fc - uj B t)]B{k, t) - aB(k, t) - ifiX(k, t). (21) 
In order to facilitate the following pcrturbative manipulation, we make the transformations: 

A{k,t) = a{k,t)e- at - i S t o dt ' [ *- Wcos{k -" Bt ' ) \ (22) 

and 

B(k,t) = b(k, t)e- at+l So «tt'[A-WcoB(fc-u, B t')]_ (23) 
Then Eqs. (17) - (19) can be written equivalently as 

^ [e at X{k,t)\ = -£[a(k,t)e- i fo dt ' [A - Wca, V'- UBt ' )] - b(k, t)e l Sl <*\±-w «»(.k-u, B t')] ]+ ^ (24) 

(Xt Z/h 

j t a(k,t) = i^X(k,t)e at+t So dt '^- Wcos ^ Bt '^ , (25) 

jb{k, t) = -ijiX(k, t)e at - l Sl dt'[A-Wco S (k-u B t')] (2g) 

As usual, we take the electrons to be initially located in band a. If the communication between the two bands is small, 
we can perform perturbative calculations by integrating Eqs. (25) and (26) and substituting the resultant expressions 
for a(k,t) and b(k,t) in terms of X(k,t) into (24), which is then integrated by iteration to the first nontrivial order. 
This gives, after summation over k, the approximate integral equation, 

p_(t) =X (t) = e- at {l + a dt'X {t')e at ' 

—/j, 2 f dt' f dt"[l + a f drX (T)e aT }J (— sin ^' f^) cosA(t' -if')) , (27) 
Jo Jo Jo \lo b 2 J ) 

where Jo is the ordinary Bessel function of order zero. Using the identity 

Jo(asinx) =J2 J 2 m {a/2)e 2mx , (28) 

m 

we can readily Laplace transform Eq. (27): 
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MP) 



l-/i 2 /(P + g) 
P + a^ 2 f(P + a) 



1 



where f(P + a) is defined by 

f{P + a)= (P + a) 2 + A 2 + 2. J ^ — ) [ 

v ' m=l 

Then (29) gives the long time behavior of Xo(i) = p-it) as 

p_(t->oo) = [l - M 2 /(a)le- t/T 



(P + a) 2 + (mule + A) 2 (P + a) 2 + (moj B - A) 2 



(29) 



(30) 



(31) 



where r, the effective decay time, is defined as r 1 = afi 2 f(a). From Eqs. (30) and (31) and the definition of effective 
scattering time r, setting tq = a^ 1 , we find 



j 2 (wy^) 

1 + (Ar ) 2 



m— 1 



1 



1 1 

- (A/Vb + m) 2 (ujBTo) 2 1 + (A/ujb -m) 2 (oj B T ) 2 



(32) 



This relation describes the leading long time dynamics of electrons. The structure of Eq. (32) shows that when A/ujb 
is an integer, t _1 /to _1 will reach a local maximum; this is illustrated in Fig. 5. These resonant peaks indicate that 
the excited state (i.e., the rungs of the upper one of the two interpenetrating Wannier-Stark Ladders) in one well is 
degenerate with the ground state of a well m unit cells away. When the scattering rate a = Tq 1 becomes smaller, 
the peaks in r 1 /^ 1 become higher and sharper (half width in A/ujb of order (ujbTo)^ 1 = a/uis), as one expects 
intuitively. Note that the decay time r in the figure is always longer than the scattering time To, even at the peaks, 
which was tacitly assumed in setting f(P + a) w f(a) at the dominant (smallest negative P) pole in (29) above. 



V. CONCLUDING REMARKS 



We have discussed the dynamics of electrons in a two-band superlattice under the influence of static electric fields, 
with Markoffian dephasing. We have given both analytic and numerical analyses of the effects of scattering and 
miniband structure on the Zener resonances and miniband localization. We found that dephasing will ultimately 
force the electrons initially located in one miniband to populate equally the two minibands. Under conditions of small 
interband coupling, and in the long time limit, we obtain the effective decay rate, whose sharp peaks are a signature 
of Zener resonances. 

Experimentally the phenomena predicted in this paper are accessible in semiconductor superlattices with wells of 
dimensions on the order of 100 Angstroms and electric fields of a few kilovolts per cm. This superlattice can be grown 
dimerizcd, with the unit cell consisting, e.g., of a pair of quantum wells, separated from the adjacent pair by a wider 
barrier than that which separates the two wells of the basis pair. Then the lowest level of an individual well is split 
within the basis pair. These two levels form a pair of minibands in the superlattice which is well separated energet ically 
from higher levels, and all but that pair can reasonably be neglected in studying the dominant electron d^namicstllE 2 ]. 
Then we can probe radiation induced by the oscillating dipole associated with Rabi flop in this systemEj. 

Work is under way to extend these results to understand the influence of scattering on the behavior of a semiconduc- 
tor superlattice in the presence of strong time periodic, as well as static, electric fields. Of particular practical interest 
are superlattices subjected to the terahertz radiation of a free electron laser. A study of the single jniniband case, 
including a discussion of chaotic dynamics in the presence of dissipation, has recently been publishedEZI by Alekseev, 
et al. 
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FIG. 1. Dependence on wave vector k of both oscillation amplitudes and the position of BO of a single component X(k,t), 
plotted as functions of dimensionless time, ust. In the upper panel (with k = 0), full Rabi oscillations can be clearly seen, 
while in the lower panel (with k — n), they are incomplete. In both figures W = (W a + M4)/2=18 meV, A = A a — A& = 20 
meV, eER — 2.09 meV, and ljb = 2.32 meV is at an avoided crossing of the interpenetrating WSL from the two bands. 

FIG. 2. The same as in Fig. 1, with W = 8.6 meV, A = 20 meV, eER = 1.84 meV, and lob has been increased to 10.2 meV 
(still at an avoided crossing). 

FIG. 3. The density matrix p~(t), with and without dissipation. We use the same superlattice parameters as those of Fig. 
1. (a) no scattering, ct/tUB = 0; (b) small scattering rate, a/u)B = 0.02; (c) large scattering rate, o/lob = 0.3. 

FIG. 4. The effect of scattering on miniband localization, (a) no scattering; (b) the complete destruction of localization by 
dephasing (with a/u) B = 0.35). Here W = (W a + W b )/2=8.6 meV, A = A a - A 6 = 20 meV, eER = 1.84 meV, and u B = 9 
meV. 



FIG. 5. Long time decay rate r 1 /r () 1 of p~(i) (see Eqs. (31) and (32)). The local maxima appear when A/cjb is an integer. 
Here fi/u B = 0.18, a/uj B = 0.2, and W/uj b = 1.3. 
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